Polariton condensation with localised excitons and propagating photons 
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We estimate the condensation temperature for microcavity polaritons, allowing for their internal 
structure. We consider polaritons formed from localised excitons in a planar microcavity, using a 
generalised Dicke model. At low densities, we find a condensation temperature T c oc p, as expected 
for a gas of structureless polaritons. However, as T c becomes of the order of the Rabi splitting, the 
structure of the polaritons becomes relevant, and the condensation temperature is that of a B.C.S.- 
like mean field theory. We also calculate the excitation spectrum, which is related to observable 
quantities such as the luminescence and absorption spectra. 

PACS numbers: 71.36.+c,42.55.Sa,03.75.Kk 
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Bose condensation of polaritons in two-dimensional 
semiconductor microcavities has been the object of con- 
tinuing experiments. Recent progress includes observa- 
tions of a nonlinear increase in the occupation of the 
ground state Q, , sub-thermal second order coherence 
|2|, changes to the momentum distribution of polari- 
tons |3jjand stimulated processes in resonantly pumped 
cavities 4j. However, unambiguous evidence for the exis- 
tence of a polariton condensate is still lacking. Theoret- 
ical predictions of the conditions required for condensa- 
tion, and of signatures of condensation, are therefore of 
great interest. We find that the momentum distribution 
of radiation from the cavity shows a dramatic change 
when condensed , which, together with the gap in the 
luminescence could provide such signatures. 

Due to the small polariton mass, a weakly interacting 
bose gas of polaritons predicts extremely high transition 
temperatures, which provides the motivation for study- 
ing polariton condensation; see [j| and references therein. 
At such high temperatures, when T c is of the order of the 
Rabi splitting, the internal structure of the polariton be- 
comes relevant, and treating polaritons as structureless 
bosons fails. This can occur at density scales much less 
than the exciton saturation (Mott) density; while current 
experiments Q are well below the Mott density, they are 
at densities where the internal structure is relevant. 

In this paper, we study polariton condensation in ther- 
mal equilibrium, including polariton internal structure. 
Although current experiments may not have reached 
equilibrium, this remains their aim; further the signa- 
tures of equilibrium condensation are instructive when 
considering the non-equilibrium case. We consider a 
model of localised excitons in a two-dimensional micro- 
cavity, building on previous studies of zero-dimensional 
microcavities 6] . This is motivated by experimental 
systems such as organic semiconductors 0, quantum 
dotsQ, and disordered quantum wells 0. Moreover it 
has been argued 01 that exciton localisation might aid 
condensation; inhomogeneous broadening due to localisa- 
tion suppresses condensation less than the homogeneous 
dephasing for propagating excitons. Finally, since the ex- 
citon mass is much larger than the photon mass, we ex- 
pect much of the physics of condensation in the localised 



exciton model to carry over to systems with propagating 
excitons, for densities much less than the Mott density. 

For simplicity, we model the localised excitons as iso- 
lated two-level systems, with the two states correspond- 
ing to the presence or absence of an exciton on a par- 
ticular site. A two-level system will be represented as 
a Fermion occupying one of two levels, further con- 
strained so that the total occupation of a site is one. 
This constraint describes a hard core repulsion between 
excitons on a single site; we neglect Coulomb interac- 
tions between different sites as being small compared to 
the photon mediated interaction. In a two dimensional 
microcavity, the photon dispersion relation for low k is 
Tiujk = Hluq + h 2 k 2 /2m, where m = (h/c)(2n/w) for a 
cavity width w. Each mode couples to all excitons, with 
different phases. Describing the coupling in the electric 
dipole gauge, the coupling strength is g — d a bV%nfiwk, 
in which we approximate ~ u>o, and d a b is the dipole 
moment. These considerations lead to the generalised 
Dicke Hamiltonian \H\ 

H = E ^ ( h ) h J - a k) + E h ^Uk 
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Here A — > oc is the quantisation area and n the areal den- 
sity of two level systems, i.e. sites where an exciton may 
exist. Without inhomogeneous broadening, the energy of 
a bound exciton is 2e = A + ujq, defining the detuning A 
between the exciton and the photon. The grand canon- 
ical ensemble, H = H — fiN, allows the calculation of 
equilibrium for a fixed total number of excitations (exci- 
tons and photons); 
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We therefore define Cbk = — A 4 an d e = e — m/2- If 
/i = 0, this Hamiltonian is an extension of that studied 
by Hepp and Lieb.0, as discussed further in 0. 

Using the characteristic density, n and energy scale 
g^/n (the Rabi splitting), the system is controlled by two 
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A(k) = AE Z + Cb% + 4SD, 



dimensionless parameters: Detuning, A* = A/g^/n and a particle in the presence of the condensate: 
Photon mass, m* = mg/H 2 ^/n. For the system in ref. Q, 
gy 7 ^ « 20meV, Tiluq « 1.6eV, and m 10~ 5 77i e . For 
delocalised but interacting excitons, assuming that the 
exciton-exciton interaction is a hard-core repulsion on 
the scale of the Bohr radius gives an upper bound on the 
density of available exciton sites as n « 10 11 riQvvcm _2 , 
where uqw is the number of quantum wells in the cav- 
ity. These values lead to dimensionless parameters m* « 
1CP 3 , and by design A* = 0. The quoted exciton den- 
sities in Ref. Q are ~ 10 9 cm -2 per well, leading us to 
estimate the fraction of occupied exciton sites to be less 
than 1CU 2 . Using these values, our work suggests (see 
Fig. that current experiments Q are in the interac- 
tion dominated mean field regime. 

Using the standard functional path integral techniques 
for the partition function as in 6] and integrating over 
the Fermion fields yields an effective action for photons: 




(6) 



At low k, the inner lines correspond to changing the phase 
of the condensate. At higher fc, unlike in B.C.S. super- 
conductivity, amplitude and phase modes get mixed. 

For small momenta, the phase mode is linear, £i = 
±cfc + 0(k 2 ). The phase velocity has the form 



c = 



J_ f 4:iM Q g 2 n\ f\ip 
2m 



N 



A po 
2m n 



(7) 




JO , 



_2_ V 



27rikr r 



_2_ V 



ipk 



We then proceed by minimising S[ip] with a static uni- 
form ip and expanding around this minimum. The mini- 
mum, ^o, satisfies the equation: 



. , 2 tanh(/3£) , 
^Wo = 9 n — Vo, 
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which describe the mean-field condensate of coupled co- 
herent photons and exciton polarisation studied in jfj . 

To consider the influence of incoherent fluctuations 
we expand the logarithm in eq. Q for fluctuations 
ij) = ipo + Sip and obtain a contribution to the action: 
±Tr(M5M- 1 M5M- 1 ). This, with the bare photon ac- 
tion, yields the inverse Greens function for photon fluc- 
tuations The resultant expression has terms 
proportional to 5 U . These arise from 2nd order poles in 
the Matsubara summation; such terms vanish on analytic 
continuation, and so do not affect the retarded Greens 
function. 

The Greens function has poles, representing the energy 
of excitations w.r.t. the chemical potential. In the normal 
state there are two poles which describe the upper and 
lower polariton dispersion 



E ± = \ ((S fc + 2~e) ± v/( 



u) k - 2e) 2 + 4g 2 ntanh(/3e) 

'(5) 

and are presented in Figure ^ The polariton disper- 
sion (JSJ from localised excitons has the same structure 
as from propagating excitons, since the photon dispersion 
dominates. In the condensed state the excitation spec- 
trum changes dramatically (see Fig ^| . There are four 
poles, reflecting the symmetry of adding or subtracting 



where the second expression is for comparison the form 
for a dilute Bose gas with contact interaction A. The 
phase velocity increases for small ipo then decreases as 
saturation of the excitons reduces the effective exciton- 
photon interaction strength. 
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FIG. 1: Condensed and uncondensed spectra, superimposed 
on the mean field phase diagram to show choice of mean field 
density and temperature. (Parameters A* = 0, m* = 0.01. 
M/SV^ — —1-4 below transition, fi/g^n — —0.6 above. En- 
ergies and densities plotted in units of g\fn and n.) 

The difference in the structure of the excitation spectra 
between condensed and uncondensed states may be seen 
in the luminescence spectra and in the angle-resolved lu- 
minescence N(k). Although four poles exist, they may 
have very different spectral weights, and would therefore 
show up differently in optical measurement. This can 
be understood more clearly after introducing inhomoge- 
neous broadening of the exciton energies, which broad- 
ens the poles. In this case, one may plot the spectral 
weight of the photon greens function, 23 [Q{iuj = z + irf)] 
as r\ — > 0, which corresponds experimentally to the ab- 
sorption coefficient. Figure shows this for the same 
parameters as figure ^ 

In the condensed case, there are three bright lines 
above the chemical potential; the Goldstone mode (which 
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is not broadened as it lies inside the gap for single parti- 
cle excitations); the edge of the gap, and the broadened 
upper polariton mode. If the inhomogeneous broadening 
is large, the gap edge peak and upper polariton will be 
broadened. This will result in a band of incoherent lu- 
minescence separated from the coherent emission at the 
chemical potential by the gap. If the broadening is small 
w.r.t. the gap then very little structure will be seen. Such 
structure can however still be seen at low densities. 
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FIG. 2: Absorption coefficient v.s. momentum (x axis) and 
energy (y axis) calculated above (top) and below (bottom) 
phase transition. Negative values indicate gain. (Parameters 
as for figure inhomogeneous broadening width 0.3g^/n.) 
The colour scale covers the range — oo to oo. 

The momentum distribution of photons may also be 
found from the photon Greens function. This is char- 
acteristic of the angular distribution of radiation; if re- 
flectivity of the mirrors were independent of angle, and 
approximating LOk ~ too then sin(#) = ck/coo- The mo- 
mentum distribution is given by: 



N(k) = lim iUl(T + f})il> h (T) 

rj— >0 \ 



(8) 



Although the system is two dimensional, a condensed 
momentum distribution can be calculated at low temper- 
atures by considering only phase fluctuations 0. Fig- 
ure |3 shows this for both condensed and uncondensed 
states. The linear phase mode leads to a power law di- 
vergence at low momenta, becoming sharper as the tem- 
perature decreases. 




FIG. 3: Momentum distribution of photons, from which fol- 
lows the angular distribution. Two uncondensed, and two 
condensed distribution are shown, for mean field densities as 
shown in the inset. (Parameters A* = 0, m* = 0.01.) 



We now estimate the phase boundary including the ef- 
fects of fluctuations. We calculate the free energy from 
the action in eq. evaluated to Gaussian order. Taking 
the total derivative of free energy w.r.t. chemical poten- 
tial gives the total density, and evaluating this at the 
critical chemical potential gives an estimate of the phase 
boundary. This follows the method of 0], but is com- 
plicated by the effects of a dynamic boson field and by 
two-dimensionality. 

The total density is 



-tanh(/?£)j 
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(dk) 2 (m) + f(Z 2 )-f(2E)+g(k)) (9) 



when condensed, with f(x) — — (nB(x) + l/2)dx/dfJ,, and 
g(k) due to 5u>; and 



P = 



1 + e 2 ^ 



K 



(dk) 2 (f(E + ) + f(E_)-f(2i)) (10) 



with f(x) = —nB{x)dx/dn when uncondensed. 

Considering the terms in eq. these can be identified 
as exciton density and coherent photon density on the 
first line. Under the integral are the occupation of upper 
and lower polaritons, and subtraction of exciton density, 
required as excitons would otherwise be counted twice. 
The cutoff K regularises the ultraviolet divergence. The 
theory is not renormalisable, due to the constant coupling 
between high energy photon modes and the excitons; but 
for low energy properties, changing the cutoff introduces 
only small corrections. 

However, as well as changing density, fluctuations 
change the critical chemical potential by depleting the 
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condensate[15|. Therefore, it is necessary to make a sep- 
arate estimate of the condensate density including fluc- 
tuations. For a two dimensional system, this can be 
achieved by calculating the superfluid response, which re- 
quires comparing the transverse (normal) response, and 
the total density: 



is just 2m. This compares to the mean field boundary, 
which at resonance is 
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where Q is the photon greens function, and p v is the 
total photon density. The transverse response may be 
safely calculated without vertex corrections, as these af- 
fect only the longitudinal response. Combining the re- 
sults of eq. © and eq. the phase transition includ- 
ing fluctuations occurs when p s = 0. As the system is two 
dimensional, one ought to consider a Kosterlitz-Thouless 
transition, at finite p s , but this introduces only small cor- 
rections. The phase boundaries for four different photon 
masses are shown in figure 0] 
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FIG. 4: Phase boundaries calculated from mean field theory, 
and including fluctuation corrections for four different values 
of photon mass. Temperature and density plotted in units of 
g\/n and n. 

At low densities, the phase transition occurs as the 
condensate is depleted by occupation of the phase mode. 
Therefore, the phase boundary has the form T c oc p/m po \ 
where m po i is the polariton mass, which, at resonance, 



T c = gVn 



2tantT 1 (I - 2p) ~ -lnp' 



(12) 



the second equality being valid for small densities. This 
temperature is effectively constant on the scale of the 
B.E.C. phase boundary. 

As the density increases, the phase boundary including 
fluctuations will approach the mean field boundary. Due 
to the slow dependence of the mean field boundary on 
density, this will occur at T ss gy/n, unless mass is very 
large. The exact temperature where this occurs depends 
on the mass. For a small mass the fluctuation bound- 
ary follows the mean field boundary over a wider range 
of densities. The rescaled photon mass, m* , provides a 
small parameter near T w gy/n, p ~ n. Physically, at 
a temperature scale of gy/n, the fluctuation modes are 
occupied to the top of the lower polariton branch, and 
occupation of the upper branch becomes relevant. 

At yet higher densities, the fluctuation result again 
differs from the mean field. For p > n the mean field 
theory will always predict a condensate. However, at 
high densities, the system is dominated by incoherent 
photons, not described in the mean field theory. In this 
limit the phase boundary is that for B.E.C. of (massive) 
photons. 

To conclude, we have calculated the spectrum of fluc- 
tuations for a model of localised excitons interacting with 
propagating photons. In the condensed state, the fluctu- 
ations about mean field include the Goldstone mode, and 
the linear dispersion causes a characteristic change to the 
momentum distribution of photons. Both the narrow- 
ing of the N(k) distribution and the gap between coher- 
ent and incoherent luminescence would provide evidence 
for polariton condensation. Fluctuations alter the form 
of the phase transition. At low densities the transition 
has the T c oc p form expected for Bose condensation of 
pre-formed Bosons. As the density increases, for small 
masses, the mean field interaction driven limit is recov- 
ered, for densities p ~ n. However, for yet higher densi- 
ties, the boundary changes again, to Bose condensation 
of weakly interacting photons. 
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